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Abstract. The approach to the consideration of the ordinary differential equations with dis- 
tributions in the classical space T>' of distributions with continuous test functions has certain 
insufficiencies: the notations are incorrect from the point of view of distribution theory, the 
right-hand side has to satisfy the restrictive conditions of equality type. In the present paper 
we consider an initial value problem for the ordinary differential equation with distributions in 
the space of distributions with dynamic test functions T', where the continuous operation of 
multiplication of distributions by discontinuous functions is defined 1171 . and show that this ap- 
proach does not have the aforementioned insufficiencies. We provide the sufficient conditions for 
viability of solutions of the ordinary differential equations with distributions (a generalization of 
the Nagumo Theorem), and show that the consideration of the distributional (impulse) controls 
in the problem of avoidance of encounters with the set (the maximal viability time problem) 
allows us to provide for the existence of solution, which may not exist for the ordinary controls. 



1. Introduction 

At the end of the nineteenth century K. Weierstrass pointed out the existence of certain problems 
in the calculus of variations where the minimum can not be achieved by any smooth function, and 
the minimizing sequence converges to a discontinuous function (e.g. the Goldschmidt solution of 
the surface of revolution area minimization problem PP). As it was further observed, the existence 
of the discontinuous solutions is not a special, but a general situation, for a large class of problems 
of optimal control theory (viewed as generalizations of the problems of calculus of variations |12|) 
with one-sided phase constraints on control |261132| . The necessity of extension of these problems in 
order to provide the existence of solutions, leads to ordinary differential equations with distributions 
or, equivalently, ordinary differential equations with measures |2l 151 12~71 1291 

Example 1 ( 20 ). The ordinary differential equation with distributions arises in the following 
optimal control theory problem 

x = x + v, x(0-)=0, x(l-) = l, v>0, (1) 

J(v) = / v(t)dt -> min, (2) 
J (0,1) 

where x is the coordinate of the space shuttle, v is the current charge of fuel, J determines the 
total charge of fuel. Problem Cpl© does not have a solution among the pairs (v,x) of locally 
summable controls v € L loc and locally absolutely-continuous x 6 AC loc , and the minimizing 
sequences converge to the elements 

v =2*0, a:(t) = | et _ 1) t>Qj (3) 
where So is the delta-function concentrated at point t = 0. 
Let us consider the following ordinary differential equation 

x = f(t,x) +g(t,x)v, (4) 

where v is a distribution, which was studied, in particular, in 0131111121112312211211, in the classical 
space V of distributions with continuous test functions jjUIJ ■ Since v S V can be represented as 
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the limit of a sequence of ordinary functions {vk}^ =1 C L loc , a solution of Q can be viewed as the 
limit 

x := lim x k in T>' , (5) 

k — >oo 

where x k e AC loc is the solution of the ordinary differential equation 

x k = f(t,x k ) + g(t,x k )v k (t). (6) 
So any solution of Q is a function of locally bounded variation, in general case discontinuous 

El EI IHI EI3 [23 EH EU 

However, the consideration of the equation @ in the space T>' leads to certain problems. 

Namely, the notations in Q are incorrect from the point of view of distribution theory, since 
Q contains the product of a discontinuous function and a distribution u 6 T>' , which is 

undefined in T>' . 

Further, let n be the dimension of the system (0}. If n > 2, then the solution JSJ) is independent 
of the choice of {w/c}^Li if and only if the condition of equality type (Frobenius condition) for 
the right-hand side of (@J is satisfied (221 - However, this condition can not be satisfied for a large 
class of problems of optimal control theory, e.g., see ^2EZ]- The approach to consideration of 
the sequence {f/c}^! as the part of the system (@J) is not new (e.g., see mEl), however, it was 
not studied from the point of view of distribution theory, since the choice of a particular sequence 
{ v k]' k= i is not motivated by any properties of the space V . 

Let us note that since V is isomorphic to the space of Borel measures (HD) on R, the problems 
described above can be reformulated for the ordinary differential equations with measures. 

The insufficiencies described above are due to the absence in T>' of the continuous operation of 
multiplication of distributions by discontinuous functions. In the present paper we overcome the 
aforementioned insufficiencies by consideration of the ordinary differential equation with distribu- 
tions in the space T 1 of distributions with dynamic test functions |17| , where the continuous 
operation of multiplication by discontinuous functions is defined (sec the discussion on the problem 
of multiplication of distributions by discontinuous functions in the classical space T>' and in other 
approached, in particular, in the Colombeau generalized functions algebra |13M14ll2*^| therein). 

Also, we study the property of viability of solutions of the ordinary differential equations with 
distributions in T . The notion of a viable solution was introduced by J. -P. Aubin |SJ[7j for the 
following ordinary differential equation 

x = f(t,x), (7) 

where / is continuous in t £ I and locally Lipschitz in x G R™. Let M C R" be a closed subset, 
to € I, fl = (ta,T) C / is an open interval, in general case unbounded. Following 0, we give the 
next definition. 

Definition 1. A solution of such that x (to) £ M and 

x(t) e M (8) 

for all t £ n is said to be viable in M (on f2). The set M is said to have the property of viability 
(on Q) for (J7J), if any solution of (J7J) such that x(to) £ M is viable in M (on ft). 

The property of viability is closely related to the problems of existence of the equilibrium points, 
construction of non-smooth Liapunov functionals, and the problems of optimal control theory 4 , 
in particular, the problem of construction of the admissible control 6 V for a controlled system 

x = f(t,x,v), x(t ) = x , veV, (9) 

such that for v = the solution of © is viable in M on £1 = (to, T*), where = max„ 6 v T (the 
problem of avoidance of encounters with the set R n \ M, or the maximal viability time problem 

EU). 

For ordinary differential equations the sufficient condition for the set M to have the viability 
property (the necessary and sufficient condition in the autonomous case) was given in |28| . 
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Definition 2 (0). Let pu '■ R™ — * R be the distance junction, pm(x) :— inf c6 M{|^ — c|}. The 
set 

K M {x) = {y G K" : Hminf (p M (z + ey)/e) = 0} (10) 
is called the contingent cone to the set M at the point x G R ra . 
Theorem 1 (Nagumo Theorem). Let M C R n be closed, M ^ 0. J/ 

f{t,x)eK M {x) (11) 

/or all t £ Q, x <E dM, where dM is the boundary of M, then M has the property of viability for 

on n. 

If the system {7J is autonomous, then Ijllfl is a necessary and sufficient condition for M to have 
the property of viability for J7J (on fi). 

After the ordinary differential equations, the conditions for viability were obtained for the differ- 
ential inclusions O E| , stochastic differential equations [5] |U| , differential equations and differential 
inclusions with aftereffect [101 123| . In the present paper we consider the property of viability for 
the systems with distributions 

x = f{t,x)+g{t,x)v, veT', (12) 

and provide the sufficient condition for a closed set M which is given by the analytic constraints 

M = {x G E n : ift(x) < 0, 1 < i < to}, (13) 

where to : M. n — > K are continuously differentiable functions (1 < i < to), to have the property of 
viability for the system (|12|l (a generalization of Nagumo's Theorem). As an application of the 
results obtained, we provide the sufficient condition for the uniform stability of the equilibrium 
points of the systems with distributions of the form 

x = f(x) + g{x)v, v G T . 

The consideration of the property of viability for the system Ijl2|) is also motivated by the 
problem of avoidance of encounters with the set K™ \ M. Namely, the problem of avoidance of 
encounters may have no solution for a system with the ordinary controls 

x = f(t, x) + g(t, x)v, v e L, (14) 

which might be viewed as the restriction of the system l|12(l to the set of regular controls v, and have 
a solution for the extended system i|12|) . i.e., when control is allowed to be distributional (impulse). 
Furthermore, the optimal solution of l|12|) has a natural interpretation as the limit of a sequence 
of solutions of l|14|) . Also, let us mention that the consideration of the property of viability for the 
system (|12|) in the space T 1 allows us to consider the trajectory at the moment of discontinuity 
(i.e., at the moment of the concentration of the delta- function at v), which is important for the 
definition of the property of viability. 

1.1. Notations. Let I = (a, b) C K be an open interval, — oo ^ a < b ^ oo. We denote by M. n 
and R nxra the spaces of the vectors and the square matrices of order n G N, respectively, with the 
elements from R, endowed with the max-norms | ■ |. Further (■, ■) stands for the inner product, (■, ■) 
stands for the componentwise product in M. n . For any p, q G R™ we have \{p, q)\ < \p\\q\- In what 
follows, we denote by <—* the embedding, i.e., the map preserving the linear and the topological 
structure. 

Let G be the algebra of bounded functions g : I — > M possessing the one-sided limits 

g(a+), g(b-), g(t+), g(t-) (15) 

for all t £ I. We identify the elements of G having the same one-sided limits, and denote the 
obtained algebra of functions (i.e., the equivalence classes of such functions) by <G = <G(7). We 
define a norm 

\\g\\ G := supmax{| 5 (t+)|, \g(t-)\}, (16) 
tei 

so G is a Banach algebra (T5|. Following ^Sji we call the elements of G the regulated functions. 
Let T(g) := {t G I : o- T (g) := g{r+) — g{r—) ^ 0} be set of points of discontinuity of a regulated 
function g G G. 
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Lemma 1 f |16|L Let g G G. Then T(g) is at most countable. 

Given g G G let us consider a partition d — {Tj}™ =1 , a < n < ■ ■ ■ < r n < b. We define 

n-l 

var/( 5 ) := sup \g(T k+1 -) - g(r k +)\ (17) 
d k=i 

(the foia? variation). If var/(<?) < oo, then we call g the function of bounded variation. We denote 
the algebra of functions of bounded variation by BV = BV(I). We define 

||0||nv:=ls(a+)|+var/(fl), (18) 

so BV is the Banach algebra. 

Let C = C(7) be the subalgebra of the continuous elements of G. Clearly, C consists of the 
bounded continuous functions on /. 

Let CBV be the subalgebra of continuous elements of BV. If g G BV, then we denote by 
g c G CBV the continuous part of g 

We denote by L = L(J) the Banach algebra of functions which are Lebesgue summable on I. 
Let us denote by AC = AC (J) the Banach algebra of functions which are absolutely continuous on 
/. By AC loc and L loc we denote the algebras of locally absolutely continuous functions and locally 
summable functions, respectively (TJJ]. We denote by BV loc and CBV loc the algebras of functions 
of locally bounded variation. 

Let us denote by F — ¥(I) the algebra of functions I — > M. We call the elements of F the 
ordinary functions. 

The definitions given above can be transferred without significant changes to the case of a finite 
closed interval. 

2. Distributions 

Let us recall the construction of the space of distributions with dynamic test functions |17j . 
2.1. Dynamic functions. Let J — [—5,5]. 

Definition 3 (J3|). A map f : I — > F(J) is called the dynamic function. 

We denote the set of dynamic functions by dW = d¥(I), and define in d¥ the pointwise operations 
of addition, multiplication by the element of R and multiplication, so d¥ forms an algebra. 

For a given / G d¥ we call /(*)(•) a dynamic value of / at t. If f(t)(-) is identically equal to a 
constant, then we call f(t)(-) an ordinary value, and denote it by f(t). 

We define the embedding F ^ d¥ as follows. For a given / G F we associate / G d¥ such that 
/(*)(•) = f(t) for alH G I. 

Let us denote 

U(f) ={t€l: /(t)(.) is an ordinary value}, D(f) = I \ U(f). 

For a given dynamic function / G d¥ and an ordinary function g G F we define the composition 
g o / G d¥ by the formula 

(9 o /)(<)(•) :=<?(/(*)(•)) (19) 
for alH G I. Accordingly, define an absolute value |/| G d¥ of / G d¥ by 

l/l(*)(0 == l/(*)(OI 
for all t G I. Let us define the support of a dynamic function / G d¥ by 

supp(/) := cl{i G 7 : t G D(/) or f{t) ± 0} G /. 

Let K C I. We define 

sup(/) := supsup(/(t)(fl)), inf(/) := inf inf ■(/(«)(«)) . 
K teK sej K teKseJ y 

A dynamic function / G cflF is said to be bounded on C /, if sup^ |/| < 00. A dynamic 
function / G d¥ is said to be non-negative on K (denote / > on K), if inf^ / > (analogously 
define a non-positive dynamic function). A constant c G R is called the right-sided limit of / at 
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t E I (denote /(r+) = c), if for any e > there exists 77 > such that sup 4g ( TT+?; ) |/ — c| < e 
(analogously define the left-sided limit at r G I = (a, b), and the one-sided limits at points a, b). A 
dynamic function / G dF is said to be continuous at r G 7, if r G U(f) and /(r) = /(r+) = / (t— ) 
(otherwise we say that / is discontinuous at r G /, denote the set of points of discontinuity of / 
by T(f) G 7; notice, that £>(/) G T(/) for any / G dF). 

Notice that if the dynamic functions in the definitions given above are the ordinary ones, then 
these definitions coincide with the ordinary definitions. 

For the purpose of construction of the space of distributions, we are mainly interested in the 
following algebras of dynamic functions. Suppose that g G dF is such that g(t)(-) G G(J) (t G I), g 
possesses one-sided limits g(a+), g(b—), g(t+), git—) for any t G /. The algebra of such dynamic 
functions (the regulated dynamic junctions) is denoted by dG, and endowed with the norm 

\\g\\ dG :=sup\g\, (20) 
1 

so dG is the Banach algebra. 

Lemma 2 f |17p. Let g G dG. Then T{g) is at most countable. 

Since D(g) G Tig), we have that any g G dG possesses the ordinary values everywhere on / 
except for certain at most countable set. For a given g G dG, we define an ordinary function g by 
g{t) := g(t) [t G U{g)). We call g the ordinary part of g G dG, and denote ord(g) := g. 

Lemma 3 ^[17]^. Let g G dG. Then g — ord(g) is defined everywhere on I for except for at most 
a countable set, is an element of G, and 

g(t+) = g(t+), g(t-) = g(t-) 

for all tel. 

Let g G dG be such that #(*)(•) G AC(J), g(t)(-l/2) = git-), g{t){\/2) = g(t+) it G I). We 
denote the algebra of such dynamic functions by sG. Then C sG dG. Let us note that for 
any g G sG we have Dig) — Tig). 

Let g G sG be such that ord(g) G BV and 

^2 v ^sej(git)is)) < 00. 

teD(g) 

We denote the Banach algebra of such dynamic functions (the dynamic functions of bounded 
variation) by sBV, and endow it with the norm 

HfllUv := \g(a+)\ + llsdlcBV + ^2 var seJ (g(t)(s)), 

teD(g) 

where g c G CBV is the continuous part of ord(.g) G BV. We have that CBV ^ sBV. We define 
analogously the algebra sBV loc of dynamic functions of locally bounded variation. 

Example 2. The Heaviside function 0@ G sBV is defined by 

<£(*) = { J; Hi] <^(r)(-)=j9(-). 
where /3 G AC (J) is such that 1/2) = 0, /3(l/2) = 1. 

Let G„, AC„, BV„, L„ and sG„, sBV rl be the spaces of vector- valued functions and dynamic 
functions, respectively, with the operations defined componentwise, and the norms agreeing with 
the norm in R™. 

2.2. Distributions. We denote by V the space of the continuous test functions, i.e., the space 
of the elements of C having compact support in /, and endowed with the standard locally-convex 
topology [301 • We denote by T the space of the elements ip G dG having compact support supp(ip) G 
/. The space T is called the space of dynamic test functions 17 . The sequence {^I^Lj is said 
to be convergent to ip in T , if (pi- — > ip in dG and there exists a compact set K C I such that 
supp(i^fe) G K for all k G N. Clearly, we have the embedding T> <— > T. 
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We use the notation V for the space of classical distributions, i.e., the space of continuous linear 
functionals V — > E [30]. We denote by T the space of distributions with dynamic test functions, 
i.e., the space of continuous linear functionals T — ► K ^7]. The value of a distribution / £ T on a 
test function ip £ T is denoted by (/, </j) £ K. The linear operations in T are defined in a standard 
way, the space T is endowed with the weak topology, so fa — > / in T' if and only if for any ip e T 
(fk,<p) -> if, <P)- 

Theorem 2 (^Tj). Any distribution in D' admits an extension from T> to T . 
Example 3. Given / £ L loc , let us define a regular distribution by the formula 

C/» := / f(W(t)dt, 



where <p £ T, <p — ord(^). Since T> <— ► T, and the linear manifold of regular distributions in 
V is isomorphic to L loc j3U] , we may identify the elements of L loc and the corresponding regular 
distributions in T 1 . 

Example 4. The delta-function 6" £ T' is defined by the formula 

(5?,<p) := J a(8)ip(T){8)ds, (21) 
where ip £ T, r £ I, and the parameter a £ L( J) such that 

a(s)ds = 1, (22) 
is called the shape of the delta-function. Notice, that for any <p £ V we have 

Or > </>) = / a(s)^(T")ds = <p(r), 



so 6" is an extension of the classical delta-function S T £ "V from 2? to T. 
For a given a satisfying (|2*2*|l we define a sequence 

<(*):=( ^W*-^))> *e(T-£,r+£), (23) 

71 w [0, otherwise, v ' 

which is called the delta- sequence having the shape a. 

We define the product of / £ T' and g £ dG by the formula 

(9f,f) ■= (f,9<P), (24) 
where gip £ T. The operation of multiplication defined by l|24|) is continuous, commutative and 
associative in the sense that (hg)f = h(gf) in T for any h, g £ dG, / £ T' |17|. 

Example 5. As follows from (|24|l . the product of the Heaviside function 9@ £ sG and the delta- 
function 6" £ T' is defined by the formula 

8*8? = U 0(s)a(s)dsj5J, 

where j(s) = a(s){3(s)/ j j (i(s)a(s)ds (s £ J) satisfies if fj P(s)a(s)ds ^ 0. 

Let / £ T 7 . Let us define the support supp(/) C I to be the minimal closed set such that for 
any ip £ T with supp((^) D supp(/) = we have (/, y>) = 0. 

A distribution / £ T' is called non-negative (non-positive) if for any ip £ T such that 93 > we 
have (/,<£>) > (or (/, ip) < 0, respectively). 

For a given distribution / £ T let us define the value of the integral over (to,ti) ([to,ti] C I) 
by the formula 



fdt:=(f,v toM ), (25) 

'(to,*i) 

where tp tQ tl (t) — 1 if t < t < t\, and <-Pt a ,U = otherwise, so <pt ,ti £ The operator of 
integration in T' which is given by (|25|l is linear and continuous |17| . 
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Let g G sBV loc , we define the derivative g G 7"' by the formula 

(g,cp):= [ 0(t)dg c (t)+ f ^(r)( S )(g(r)(s)) s d S . (26) 

where tp G T, g c G CBV loc , and the set of points of discontinuity T(g) is at most countable by 
Lemma|21 If for any r G T(g) a>(<7) = <?(t+) — s(t— ) ^ 0, then 

<7 = <7c + E "rla)*? ( 2? ) 

in 7"', where the shape of a delta-function a T (s) = (g(r)(s)) s /a T (g) (s G J). 
Example 6. The derivative of the Heaviside function 0£ is the delta- function 6", 

qP _ S a 

where a = [3. 

In what follows the notations V' n and 7^ stand for the spaces of the vector- valued distributions 
with the convergence and the operations defined componentwise. 

3. Systems with distributions 



Let D C I x R™ be open. Consider in 7^ = 7^(7) an initial value problem 

x = f(t, x) + g(t, x)v, x(t -) = x , (28) 

where (to, Xq) G D, 

1) The function /:£)—> W 1 is Caratheodory in (t, x) G D and Lipschitz in a; with the constant 
K f > in D. 

2) The function g : D — ► R™ xn belongs to sG„ in i (in particular, <? is continuous in i) and 
Lipschitz in x with the constant K g > in 7). 

3) The distribution v G 7^ is defined by 

v = ii, 

where it G sBV^ 30 , we assume that for any t G T(u) <j t (u) ^= 0, so according to the definition of 
the derivative in 7^ 

v = ii c + (vt(u),6^), (29) 

where S" T G %[ is the vector- valued delta- function, (8^ T ,tp) := ((&t t , (fj) -_ v where tp G 7", 
a T = (c4_) " =1 G L n (J), (, ) is the componentwise product in W 1 . 

Let us note that g may be discontinuous in t. This is important, since in the problems of optimal 
control g may also depend on the (discontinuous) ordinary control (in t), which is not considered 
here. 

A solution of the initial value problem l|28|) on £1 C I is the dynamic function x G sBV[° c (Sl) 
such that x(t)(s) G D (t G fi, s G J) and ggj is satisfied in 7^(fi). An ordinary part £ G BVj° c (ft) 
of the solution x is called the ordinary solution of i|28|) . Let us note that in contrast to the 
approach based on the space of distributions T>' the operations of multiplication, differentiation 
and composition arising in l|28|) are correctly defined in the sense of the distribution theory in 
7^(f2) (see the Introduction). 

In the formulation of the next theorem we put Q = I. 

Theorem 3. Let x G sBV 1 " be the solution of (SgJ, x G BV 1 " &e i/ie ordinary solution of 
TTien 

= x + / /(r, x(r))dr + / <?(r, i(r))c?u c (r)+ 
Jto J t 

^( 7 r(l/2) - 4(r-)) - E (7r(l/2) - i(r-)), (30) 

r<t r<t 



8 



D. KINZEBULATOV 



and the dynamic value 7t (-) := x(t)(-) satisfies 

7t(s) = g(T,<y r (sj)(s)(a T (u),a T (s)), j T (-l/2) = x(t-), (31) 

where T(u) = {r} is at most countable by Lemma\Q g(t,x)(s) is a dynamic value of g at (t,x), 
x(t+) = x(t+) = 1t (1/2). 

Conversely, any x £ sBV^ satisfying <|3U|) . i|31[l is a solution of (|28p. 

Proof. Observe that by our assumptions a solution of the problem l|31fl exists and is unique. Let 
us show that the second statement holds. Suppose that there exists a solution x £ BVj^ c of the 
integral equation H3()(l . Then we may define x £ sBV loc as in the formulation of the theorem. First, 
show that x £ sMY° c satisfies the initial condition in H28JI. We have that x(to— ) = x(to— ), the 
limit x(to— ) exists since x £ BV„. Since all the integrals in 130(1 are continuous in t, we have 



lim / f(r,x(r))dr = 0, lim / g(r, x(r))du c (r) = 0, 
t->to- J tg t^t - J to 

so the limit of the step part in (|30[1 is 

Jim £( 7t (1/2)-s(t-)) - ^ (7.(1/2) -f(r-)), 



t->t - 

T<t T<t 



i.e., x(t — ) = xq. Second, we show that x £ sBV^ c satisfies the differential equation in (|2"%|l in 7^. 
According to the definition of the derivative, we have 

(x,<p) = I <p(t)dx c (t) + I V>(T)(s)«r)(s)) s ds, (32) 

where x(t)(s) = 7 r (s), so differential equation l|28(l is equivalent to 



U(t)dx c {t)+ / vW(s)Wt)( s )); 



ds 



tGT(u) 



f(t,x(t))tp(t)dt + g(t,x(t))du c (t)+ / f(T)(s)g(T,x(T)(s))(s)(a T (u),a T (s))ds. 

Let us show that the equality 



ip(T)(s)(x(T)(s)) s ds = / ip(T)(s)g(T,x{T)(s))(s)(a T (u),a T {s))ds (33) 
J J 

holds for any <p £ T, r £ T(u), and 

± c = f{t,x) + g(t,x)ii c . (34) 

Indeed, the equality (fTTT^f) follows from lf3"T)> . where x(r)(s) = 7 T (s) (t £ T(uj). Further, by the 
definition of the regular distribution in 7^, the equality (|34|l is equivalent to the equality 

tp(t)dx c (t) = £ <p(t)f(t, x{t))dt + J <p(t)g(t, x(t))du c {t) 

for any ip £ T n . Then the last equality is equivalent to 

x c {t) = x c (t ) + / f(r,x(r))dr+ f g(r,x(r))du c (r) (35) 



for all t £ I. As follows from the construction of x £ sBV loc in the statement of the theorem, the 
value of the jump of the solution x(t+) — x(t—) can be found from the dynamic value x(t)(-) = 
7t-(-). Consequently, the equality of the continuous parts (|35|) implies (|3U|l . Since the equalities 
obtained are equivalent to the equalities (|33[L((34|) in T„, the dynamic function x £ sMY l ° c is a 
solution of the initial value problem l|28l) . 

We now show that the first statement holds. Let x £ sBVj° c be the solution of the problem 
(J2BJ. Then (JSSJ and are true. Further, the equality (E3J is equivalent to 

<p(t)dx c (t)= ( 0(t)f(t,x(t))dt+ { 0(t)g(t,x(t))du c (t) 
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for all <p £ T n . According to DuBois-Reymond Lemma we have 

x c (t) = x c (to) + [ f(r,x(r))dr+ f g(r, x(r))du c (r) (36) 

for all t £ I. Analogously, as follows from DuBois-Reymond Lemma, the equality (|33[l implies that 
(|31[) holds for any r £ T{u). Then (|30fl is true, as follows from the definition of the space sBVj° c 
and the initial condition x(t$— ) = xq. □ 

Theorem 4. There is a constant h > such that there exists a solution x £ sWV n (to —h,to + h) of 
the problem ()28[1 , which is unique in the sense that it coincides with any other solution of (|28|) on 
the common interval of definition in (to — h, to + h), and depends continuously on u|( to _/ l to+ / l ) £ 
sMY n (t -h,t Q + h). 

Proof. Let us denote Ih ■— (to — h,t n + h), if := (to, to + h), 1^ := (t — h, to) C I = (a,b). Let 
us consider in W n (Jt") the following integral equation 

x(t)=x + [ f(r,x(r))dr+ f g(r,x(r))du c (r) + ^ ( 7r (l/2) -4(t-)), (37) 

^0 ^0 t <T<t 

7T(a)=ff(r,7r(s))(*)<ov(«) J a T (*)>, Tr(-l/2) = x (r-), (38) 

where r £ T(u) C it , the dynamic function u is assumed to be restricted to if . Clearly, if the 
solution x £ MV n (lt) of (|37|l l|38|l exists, then x (to+) — Xo. Let us show the existence of solution 
of in BV„(/+) for some h > 0. Let N > 0. We define 

M' g l = max{|fl(t,a;)(s)| : s £ J, t £ [t ,t + h], \x - x \ < N} > 0, (39) 



Q+ = {x£ BV„(/+) : \\x - a;o|lBv n (j+) < (40) 

so that [to, to + h] x {x £ R" : \x — x \ < N} C D. Obviously, is a complete metric subspace. 
We define a map P on by the formula 

P(x)(t)=x + f f(r,x(r))dr+ f g(r,x(r))du c (r) + ]T ( 7t (1/2) - £(r-)), 

J *<> Jto t <T<t 

where t £ if, x £ Q\. Then 

\\P(i) - ^o|| BVn (/+) < / l/MM) \dr + M^(var +( w )) 2 + M ff ft var J+ (u). 

h 

Then there exists h > such that ||P(x) — 2 ; o|Ibv„(/+) ^ ^> i- e -' ^ : Qfr ~^ ^fe- Let us snow that 
there exists h > sufficiently small such that there is A > such that 

\\P@) " P(v)\\mv n (i+) < M\* - y\\mv n (i+) ( 41 ) 

for all x, y £ Q^. Indeed, to have the inequality (|4*T|) satisfied it suffices to put 

A = Kfh + K g va,Tj+ (u)h + K g vaXj+ (w)var* o (u) > 0, (42) 

so we can find h > sufficiently small such that A < 1. Then according to the Fixed Point 
Theorem |25| the mapping P has the only fixed point in Qt, which is the only solution of the 
integral equation (|37[H|38[1 . Since the solution -f T of the problem (|3I() exists on J for any r £ T(u), 
and T(u) = T(x) = {t £ if : x(t+) ^ x(t-)}, we may define x £ sWV n (I h ) by x(t) = x(t) if 
t T(u), and x(t)(s) = 7t(s) (s £ J) if t £ T(u). According to Theorem|21 the dynamic function 
x £ sWV n (l£) is the solution of the initial value problem (|28|l on if. 
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Further, by the change of independent variable in (|37[) from t to — i we can show that there 
exists h > such that [to — h, to + h] x {x G E™ : \x — xo\ < N} C D and the following integral 
equation 

x(t)=x + / f(r,x(r))dr + / g(r, x(r))du c (r)+ 

(lr(l/2)-x(r-))-YhAl/2)-x(r-)), (43) 

t -h<T<t T<t 

endowed with condition (|38[l . where t G T(u) C /h, and the dynamic function ti is assumed to be 
restricted to 7^,, has a unique solution in 

Q h = {xe BV„(4) : p - zolU„(/ h ) < W}. 

Now let us show that the mapping P = P u depends continuously on u 6 sRV n (l£). Suppose 
that we have {u k }f =1 C sMV n (l£), u G sBV„(/,|) and 

Ufe — > u 

in sBV ra (7 ? |"). Let us show that for any x G we have P Uk (x) — > P u (x) in Q^. Indeed, we have 

||-P«(£) - f« fc (^)||]gv(/+) ^ M^vaxj+(u - u k )+ 

M var /+ (w — u fc )var /+ (w) + M^var /+ (M fc )var /+ (u - u k ), 

" h h " h h 

where {var / + (wfc)}^! is bounded, so 

\\Pu{x) --P« fc (^)|| B v„(z+) ~^0as ||w-Ufc|| sB v„(i+) 

As is shown in |25| . the inequality above implies the convergence x k — > x in Q^, where G is 
a fixed point of the mapping P Uk (k G N). Consequently, the solution x G BV„(i^~) of the integral 
equation H3U|) depends continuously on u G sBV rl (/,+). 

Let us show that the solution x G sBV,^/^) of the problem (|28Jl depends continuously on 
u G sBV„(/ ft h ). Let us denote the dynamic values (i T (-) = u(r)(-), /?£(•) = Ufe(r)(-) G AC(J) 
(t G and denote the solution of the system (|31|l corresponding to u k by 7.™ G AC(J). 

Employing the known estimations we have 



varj( 7r - 7*) < M s ft l |o>(u) - CT T (u fc )|var, 7 (/3 T ) + |a T (u fc )|varj(/3 T - + 

lf 9 |o- T (ufc)|e n ^ iK « h Af*fvar J+ (u)|o- T (u) - u T {u k )\ + |a T K)|var /+ K)varj(/3 T - $) 

^ \ h h 

where the sequences {var / +(itfc)}, {|cr r (ttfc)|} are bounded. Consequently, 

h 

IN - X \\aBV n {I+) ~* if IK ~* U ILbV„(/+) ~» °- 

The case of the left-sided neighborhood 17 is treated similarly, so the solution x G sMV n (Ih) of 
the initial value problem l|28|) depends continuously on u G sWV n (Ih)- D 

The existence of a non-continuablc in D solution of the initial value problem (|28[) can be shown 
in a standard way. 

Let D = I x N, N C K n be open. An initial value problem 

x = f(t, x) +g(t,x)v, x(t -)=%o (44) 

where v = u G V n , u G BV^ oc , was considered, in particular, in E01 E3 HD EU The solution 
a: G BVj°° of the initial value problem (|44|) is defined by 

x := lim x k in V' nl 

k — >00 
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where x k G AC loc is the solution of the initial value problem for the differential equation with the 
ordinary right-hand side, 

x k = f(t,x k ) + g(t,x k )v k (t), X k (t ) = Xq, 

where v k = u k G Lj° c , u k G AC^ C , and v k — > v in V' n |271 I29| . If the function g is continuously 
differentiable in (t, x) G .D, then the necessary and sufficient condition for the uniqueness of the 
solution of i.e., its independence on the choice of the approximating sequence {wfc}^L 1; is the 
Frobenius condition 

[g m ,g l ] x (t,x) = 0, (t,x) eD, l<l,m<n, (45) 

where g m is an m-th column of g, and [•, -] x is the Lie bracket in the variable x, e.g., see [23 - If 
the condition (|45|l is satisfied, then the solution x G BV|^ C of the initial value problem (|44ll satisfies 



x(t) = x + / f(r,x(r))dr + / g(r, x(r))du c (r)+ 

Jt Jt 

^2( 7t (1/2)-x(t-)) - £( 7t (1/ 2 ) -*(,--)), (46) 

T<t T<t 

7r(a)=s(r,7r(a))o-r(«) ) 7r(-l/2) = z(t-), (47) 

where T(u) = {r}, x(r+) = 7(1/2) 

As is well known, the condition l|45|) is a necessary and sufficient condition for the invariance of 
the value 7(1/2) G M. n with respect to the choice of the function a G L„(J) satisfying (|22l) . where 
7 is a solution of the problem 

7(«)=$0r,7(*))a(«), 7(-l/2)=7o, (48) 

where r G /, 70 G N [22] - Consequently, H45|) is equivalent to the condition of the independence of 
the ordinary solution of problem (|28l) on the choice of the shapes of delta-functions a T (r G T(u)) 
(see the Introduction). 

Let us note that if the condition l|45|l is satisfied, then l|47|l coincides with (|31|l for a k = 
(1,...,1) T - 

Example 7. Let I = (— 1, 1), D = I x M. Let us consider in T' the initial value problem, 

i = ^a:, x(-l/2-) = x 0) (49) 
where xq G R. According to Theorem|31the solution x G sBV of the problem l|49|) is given by 

»(*) = {Z, Itl m(s)= X oe^^ (seJ), 
where J = [-1/2,1/2]. 

Note that the same value of the jump of the solution at t = can be obtained if the delta- 
function Sq in l|49(l is replaced by the terms of the delta-sequence l|23|l having the shape a. Also 
note that the value of the jump of the solution at t = is independent of the choice of the shape 
a since for n = 1 the Frobenius condition H45|) is always satisfied. 

4. Viability Theorem and stability analysis 

Let f2 = (to, T) C /, let M C W 1 be a closed subset, D = I X N, where N is an open subset, 
M C N. The following extends Definition ^ 

Definition 4. A solution of the system (|28|l such that x(to—) G M and 

x(t)(s) G M 

for allt<EQU {to}, s G J, where to G /, is said to be viable in M on f2. The set M is said to have 
the property of viability for (|28|l on Q, if any solution of (|28|l such that x(to— ) G M is viable in M 
(on SI). 



12 



D. KINZEBULATOV 



Suppose that the set M C R" is given by 

M = {x e K" : rn(x) < 0, 1 < i < m}, (50) 

where r\i : M™ — > K. are continuously differentiable on M™ (1 < i < m), so M is closed in M". 
Clearly, for each l<z<mifa:S R™ is such that r\i(x) = 0, i)i(x) ^ 0, then 

: m (p)<o}(x) = {y : (vi(x),y) < 0}. 

Then according to p. 224] for any x G <9M such that ?)i(a;) el" (z G L x = {i : rn(x) = 0}) are 
linearly independent, the contingent cone Km(x) is given by 

K M {x) = {f/eK": {m{x),y) < 0, i G L x }. (51) 

The following theorem follows immediately from the Nagumo Theorem. 

Theorem 5. Let f]i(x) G R™ (i G L x ) be linearly independent for any x G dM. If 

(rn(x),f(t,x)) <0 (tGia) 

/or all t £ fl, x £ dM , then M has the the property of viability for J7J) (on Q). 

Let us consider in 7^ the following differential equation with distributions of the form (J2SJ, 

oo 

i = /(t,a:)+5(t J a;)t; 1 « = «; + C> e T n, (52) 

fc=i 

where the function w : J — > R™ is continuous, {ife}^^ C 7, <5" fc fc G 7^, the coefficients G R™ 
are such that the primitive of v in T 7 [ is in sBV[° c , G C„(J), the functions /, g satisfy 1), 2) 
and / is also continuous in £ 6 S7. We suppose that for any xq G D there exists a non-continuable 
solution of the initial value problem for (|52|) in D with the initial value x(to — ) = xo- 
The following statement generalizes Theorem 

Theorem 6. Let f)i(x) G R™ (i £ L x ) be linearly independent for any x G <9M. J/ 

(r)i{x),f(t, x) + g(t, x)v) < (i G L x ) (53) 

m 7"'(f2) /or aZ/ x G <9M ; i/iera M has the property of viability for 152|) (on Q). 

The definition of a non-positive distribution in T was given in Section 2. As follows from 
(|53|l and the examples below, the property of viability depends on the choice of the shapes of the 
delta-function in v, including the case where the Frobcnius condition (|45|l is satisfied. 

Proof. 1) Let x G dM, i G L x . Let us show that if the inequality (|53[1 holds, then 

(f h (x),f(t,x)+g(t,x)w(t)) <0 (54) 

for alH G O, and 

(%(£), S( r fc,a0(.s)(cfc,afc(s))) <0 (55) 

for all s G J, A; G N. According to the definition of a non-positive distribution in 7"', the inequality 
(|53|l implies that 



(m(x)J(t, x) + g(t, x)w(t))tp(t)dt + YJfH, g(T k ,x)(s){c k , (*°* , if))) < 
r fc=l 

for any (/? G T, > 0, where tp — ord(^), i.e., by the definition of the delta- function 

OO , p \ 

(i)i(x), f(t,x) + g(t,x)w(t))ip(t)dt + (fji, / g(r k ,x)(s) (c k ,a k (s)(p(T k )(s)) ds \ <0 (56) 



fe=i 



for any tp G 7", </? > 0. Let fen £ N be fixed. Let (/? G T be such that tp = 0, ¥'( T fc)(') = (fc ^ fco)> 
vC^oX*) ^ 0- Then (jHEJ implies that 



/ 9( T k B ,x)(s) {c ko ,a ko (s) tp(r ko )(s)) ds ) < 0. (57) 
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Due to the linearity of the integral, the inequality l|57|l is equivalent to 

(r]i(x), g(r ko , x)(s)(c ko , a ko (s)))f(T ko )(s)ds < (58) 

for any <f{T ka )(-) > 0. Then due to the continuity of the functions in l|58[l we have 

(fii(x),g(Tk ,x)(s)(c ko ,a ko (s))) < 

for all s G J. Since ko G N was chosen arbitrarily, we obtain l|55|) . 

Let ip G T, ip > 0. Since the change of the dynamic values of if in finitely many points T k does 
not change if = ord(<,5), and the series in the right-hand side of l|56[) converges, we obtain that 

J (fii(x), f(t, x) + g(t, x)w(t)) ip(t)dt < 

for all if G T, if > 0, which implies 1)54(1. 

Since x G dM, i G L x were chosen arbitrarily, according to Theorem [S] we obtain that the 
conditions l|54|l and (|55|l imply that M has the property of viability for 

x = f(t,x)+g(t,x)w(t), (59) 

and for the system i|31|) (for any k G N), respectively. 

2) Consider first the particular case where there exists / > such that r k+ \ — r k > I for all 
ken. Let ft = (0,oo). 

Let Ti > 0. Then x(0) = Xo and by Theorem [21 since Tfc+x — r k > I > for all fc 6 N there exists 
rj > such that a; has the ordinary values x(t) and 

x(t) = f(t, x(t)) + g(t, x(t))w(t), (60) 

for all t G [0, rj). As is mentioned above, M has the property of viability l(60|l on (0, 77), so cc(t) G M 
for all * G [0,77). 

Now suppose that t\ = 0. Then by Theorem |3| the jump of the solution x at t\ ~ can be 
found from (|31|l at point T\ . Due to the remark above we have that M has the property of viability 
for at the point t\, so since 7i(— 1/2) = x(0— ) G M we have that 7i(s) G M for all s G J. 
Consequently, a;(0+) = 7(1/2) G M, Analogously to the first case we obtain that there exists 
7] > such that x(t)(s) G M for all t G [0, 77), s G J. 

We show that x(rj)(s) G M for all t > 0, s G J. By the change of the independent variable t we 
obtain that the inclusion x{t^ — ) G M implies that there exists r\ = ri(to,x(to—)) > to such that 

x{t)(s) G M (61) 

for all t G [to,i]), s G J. Thus, we obtain a strongly monotonically increasing sequence {t k } k x L 1 
such that M has the property of viability on (0, t k ) Suppose that the sequence {ifcjfcLi is bounded 
from above. Consequently, t k — > t* from the left, where < t* < 00. By our assumption the 
solution x is defined for all t > 0. Since M is closed, the limit x(t*—) G M. Thus, we may change 
the independent variable t, and apply the same arguments for to = t* . As a result, we obtain a 
contradiction with the assumption that t* is the maximal possible, so (|61|l holds for all t G (0, 00), 
s G J. 

The case of bounded f2 is treated similarly. 

3) Consider the general case. Without loss of generality we give a proof for the case f2 = (0, 00). 
Let [c, d] C fi, u e sBV„(c, d), v = ii G 7^(c, d), 

OO 

k=l 

where q : I — > R™ is continuously differentiable, w — q, p k £ AC„(J), Qfe = (fc G N), and 
without loss of generality G (c, rf) (fc G N). Let us define 

\c h \>yi 
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where vi = ui G T^(c,d) contains the linear combination of delta- functions. Since 



\\ u - u l\\sMV n {c,d) 



-> 



|CA,|<1/1 

(I — > oo), by Theorem0]we have the convergence x; — > a; in sBV n (c, d) of the sequence of solutions 
xi G sBV„(c, d) of the initial value problems (|52(l for = u;. We apply the results obtained above 
to the initial value problems (|52[1 corresponding to vi G 7^, so, any solution x; is viable in M on 
(c, d), i.e., xi(t)(s) G M for all t £ (c,d), s £ J. Convergence x; — > x in sBV„(c, d) implies that 

xi(t)(s) -+x(t)(s) (t G (c,d), s G J). 

Since M is closed, we have x(i)(s) G M for all i G (c,d), s G J. Since [c, d] C f2 was chosen 
arbitrarily, we obtain that x(t)(s) G M for all t > to, s <E J. □ 

Corollary 1. Lef ?7i(x) G l n (i G L x ) be linearly independent for any x G dM. If 

(f H (x),f(t,x)+g(t,x)w(t)) <0 (62) 

/or all t x £ dM, 

(fii(x),g(T k ,x)(s)(c k ,a k (s))) <0 (63) 
/or all s £ J, k E N, x G 9M, t/ien M /ias t/ie property of viability for (|52|) (on O). 

Proof. The proof follows from the proof of Theorem □ 

Example 8. Let I = (-1, oo), D = I x (-1, 2). Let O = (0, oo) C /, {ife}^ C fi, r fc -> oo. Let 
us consider in T the following ordinary differential equation with distributions of the form Q52p. 



2 x y ^ " rfc 
7 fe=i 

Let r)(x) = {x- 1/2) 2 - 1/4, so 

M = {x 6l: ?7(x) < 0} = [0,1]. 

Notice, that ?)(x) G {— 1, 1} if x G <9M = {0, 1}. Since the system (j^U) has form (gU for /(x) = -x, 
,g(x) = 1/2 — x, and v = Y^kLi ^ > we ma y a PPly Corollary^ We have 

< 0, »)(Q)/(0) < 0, 

and 

fi(l)g(l)a k (s) < 0, 7j(0) 5 (0)a fe (s) <0 

for all * G J, fc G N, where t)(1) = 1, g(l) = -1/2, t)(0) = -1, g(l) = -1/2, /(l) = -1, /(0) = 
and afc > (s G J, fc G N). According to Corollary ^ the set M has the property of viability for 
the system (I64f) on ft. 

Let us consider the applications of Theorem to stability analysis. 

Let I = (a, oo). Consider the ordinary differential equation with distributions of the form l|28|l. 

x = f(x)+ g(x)v, (65) 

where / : D — > W l , g : D — > K™ xn are Lipschitz in £>. The solution of l|65|> which is identically 
equal to a constant x* G D is called the equilibrium point (clearly, x* is an equilibrium point of 
J2U if and only if f(x*) = g(x*) = 0). 

Following the standard terminology, we note that the equilibrium point x* G D is uniformly 
stable, if there exists a sequence of the closed subsets M n+ i C M k C D (k G N), where M k 
contains an open neighbourhood of x* G D and is also contained in another open neighbourhood 
of x* G D, such that for any I G N there exists k ^ I having the property: any solution x of (|65|l 
with x(to— ) G M k is viable in Mi on f2 = (to, oo) c 7. 

Theorem 7. Suppose that x* G 7) is an equilibrium point of Q65[). and 

(x-x*,/(x)+5(x)w) <0 (66) 

in T^ifi) for any x G D such that \x — x*| 2 = 1/i, where \ ■ \ 2 is a Euclidean norm in K™ , Z G N. 
TTien x* is uniformly stable. 



<e, a fc >0, « fc £C(J). (64) 
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Proof. Let us note that in the definition of the uniform stability above we may have I = k, that 
is, it suffices to prove that for each I G N the subset Mi has the property of viability for (??) on 
Q = (to, oo). Let 

Mi = {x e D : Mx) < 0}, where (5Ax) = \x -x*\\- i 

V 

Then the sequence {Mi}f^ 1 satisfies the conditions above. We have 

fax) = 2(x - x*) ? 

for any x € D such that \x — a;* 1 2 = l/l. According to Theorem the inequality (|<jfc>|) implies that 
Mi possesses the property of viability for (|65|) on fi, so x* is uniformly stable. □ 

Analogously to Theorem we obtain the following corollary of Theorem 

Corollary 2. Suppose that x* € D is an equilibrium point of (|65l) . and 

(x-x*,f(x)+g(x)w(t))^0 (67) 

for all t G fi = (iooo) and all x & D such that \x — x*\2 — l/l, 

(z - x*,g(T k ,x)(s)(ck,a k (s))) ^0 (68) 

/or a^/ s £ J, t 6 fj and all x & D such that \x — x*\2 = l/l, where | • |2 is a Euclidean norm in 
R n , I G N. Then x* is uniformly stable. 

Proof. The proof follows from the proof of Theorem □ 
Example 9. Let us consider the ordinary differential equation with distributions 



1 \ 00 

o- EC (69) 



fe=i 

where / = (—1, 00), the shape a k is continuous. Clearly, l|69|l has an equilibrium point 

1 

** = -. 

The equation l|69[) has the form of l|65|l for f(x) = 0, g(x) = \ — x, w(t) = 0. We have to show 
that the inequality (|68|l holds. Indeed, the inequality (|68|) is equivalent to 

^ - x\ a k (s) (x - ^ = - Q - x^j a k (s) < 

for any s € J and any x such that \x— 1/2| = l/l (I G N). According to Theorem0the equilibrium 
point a;* = 4 of equation l|69|l is uniformly stable. 

5. The statement of the impulse problem of avoidance of encounters 

Let M c M. n be a closed subset, D = I x N, where N is an open subset, M C N. Let us 
consider in K™ the following controlled system of the form (0 . 

x = f(t, x) + g(t, x)v, x(ta) — xq, v g V, (70) 

where /, g satisfy conditions 1), 2), the function / is continuous in t, and the set of admissible 
ordinary controls V is given by 

V = (v G L„(7) : > 0, Ui(s)ds < V, 1 < * < n| , (71) 

where = (i 0) 7 1 ) C I, V > is given, v = (v;)?=i- 

Let Vm(? 1 ) G V be the maximal set of admissible controls such that 

v G Vm{T) implies that x is viable in M on O = (to,T) C 7, 

where x G AC ra is the solution of system l|70|) . According to |U > we cau the following maximization 
problem 

T->max, v£V M (T) (72) 
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the problem of avoidance of encounters with the set R™ \ M. 

Along with the system J7UJ), let us consider in 7^ the following controlled system 

x = f(t, x) + g(t, x)v, x(t — ) = x Q , v S V', (73) 

where /, g satisfy 1), 2), the function / is continuous in t, and the set of admissible distributional 
(i.e., impulse) controls V" is given by 

V = ju e : v, > v r ds < V,l < i < nj , (74) 

where V > is given, v — (vi)™ =1 (the definitions of a non-negative distribution and integral of a 
distribution were given in Section 2). Notice that V C V, where the elements of V are viewed as 
the regular distributions in T 7 [. Thus, system (|75|) (|7l*|l is an extension of system {7UJ)([72). 
We define V' M (T) C V" to be the maximal set of admissible controls such that 

v G V' M (T) implies that x is viable in M on f2 = (to,T) C I, 

where x £ sWV n is a solution of system (|73f) . Analogously to [21] , we call the maximization problem 

T^max, veV' M (T) (75) 

V 

the impulse problem of avoidance of encounters with the set M" \ M. 

As follows from the next example, the problem of avoidance of encounters may have no solution 
(T*,v*) for the system (I7U)) (1711 . but may have a solution for the extended system (|7H1) (|7^|) . 

Example 10. Let / = (—1, oo), D = I x R. Let r/(x) = x 2 - 1, so 

M = {x € M : n{x) < 0} = [-1, 1]. 

Let us consider in T the following controlled system 

x = x-v, »eV, x(0-) = l, (76) 

V = jv E T : v > 0, J vdt<^ 

Let us show that the solution of the problem of avoidance of encounters with the set R \ M for 
(123) is given by £1* = (0, T*) C 7, 

T*=ln(2), v* = ^, a>0. 

Suppose that v = 1/25q, a > 0. Then according to TheoremO^we have x(0)(s) = 7(5) (s 6 J), 
x(t) = l/2e', t e (0,ln(2)), where 



< 7 (s) =1-1 J ^ a( V )d V < 1 



for all s G J, so x is viable in M on SI = (0, ln(2)) C I. Along with that, for any regular control v 
solution of H76|) is given by 



e 



x(t) = e* - ^- / w(s)e _s ds 



for all £ > 0. As follows from the obtained representation, there exists £ > such that x(t) > 



z t_ e te 



-t 



for alH > 0. Thus, 



2 

T<ln (jzW) <M2), 

so T* = ln(2) is the maximal viability time. 

In the subsequent paper we provide the necessary conditions for optimality in the impulse 
problem of avoidance of encounters with the set R" \ M. 
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